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Abstract. The cross-slip process of a screw <a> dislocation from the basal to the 
prismatic plane in magnesium was studied using the density functional theory and 
the molecular dynamics calculations. An atomistic method for calculating the total 
Peierls energy map has been devised to track the transition path of a dissociated 
and/or constricted screw <a> dislocation in the cross-slip process. The barrier of a 
screw <a> dislocation from the basal to the prismatic plane is estimated by the density 
functional theory for the first time to be 61.4 ± 2.0 meV per Burgers vector length. 
The activation enthalpy for the cross slip is calculated using a line tension model 
based on the density functional theory to be 1.4 to 1.7 eV, which is in reasonable 
agreement with experiments. On the basis of the results, the effect of temperature 
on the cross-slip process of the dissociated screw <a> dislocation on the basal plane 
is studied in detail using the molecular dynamics method with the embedded-atom- 
method (EAM) interatomic potential, in which the critical resolved shear stress for 
the cross slip is evaluated. It is confirmed that the bowed-out dislocation line on the 
prismatic plane consists of slightly dissociated rectilinear segments with connecting 
jogs at low temperatures and, as the temperature rises, the curved dislocation line 
becomes smooth with many segments. The motion of an <a> dislocation on the 
prismatic plane is jerky in the low temperature region, while it is retarded by the 
formation of the largely dissociated plateau segment above the room temperature. A 
large reduction of the critical shear stress for the cross slip is obtained when the <a> 
screw dislocation interacts with a hard-sphere particle placed on the basal plane in the 
low temperature region. 


1. Introduction 

Magnesium alloys are the key candidates for the development of lightweight structural 
materials. Formability of these strongly anisotropic hexagonal metals at room 
temperatures is the stumbling block to achieve this goal, which is to enhance plasticity, 
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in particular, near the <c> axis direction to induce a general homogeneous deformation 

Cl¬ 
in pure magnesium, the primary slip system is basal and the secondary systems 
are prismatic and pyramidal along with twinning. The <a> and <c+a> are the two 
Burgers vectors which induce non-basal deformations[2]. Among them, there are two 
slip systems that are observed and studied: one is the prismatic slip of <a> screw 
dislocations by the cross slip from the basal plane, and the other is the pyramidal 
slip of <c+a> dislocations. The mechanism of the formation of a <c+a> dislocation 
from the combination of <a> and <c> dislocations through two cross-slip processes 
is proposed[3), but not fully confirmed. It is known that ductility of magnesium is 
dramatically improved by adding the elements of Li, Y, Ce, etc. 0 U 0 U\. According 
to these transmission electron microscopy (TEM) observations after deformation, a 
high activity of <c+a> dislocations is observed in Mg-Li and Mg-Y alloys, while 
the prismatic <a> slip is activated in Mg-Ce alloy. The results depend on the solute 
content and the range of temperature. The cross slip of <a> dislocations to non-basal 
planes is also observed in deformed Mg-3Al-lZn-0.2Mn(AZ31B)[8j|. Thus the prismatic 
<a> slip and the pyramidal <c+a> slip are two important processes for the non-basal 
deformations. Here we focus on exploring the mechanism for the cross-slip process of a 
screw <a> dislocation from the basal to the prismatic plane in pure magnesium. 

The tensile deformation experiments are performed on single crystals of magnesium 
and Mg-Li alloys to study the effects of temperature and solute concentration on 
prismatic slip [9], |10[ [Til I2j- In particular, the in situ experiment of prismatic glide 
in magnesium single crystal is conducted by inhibiting the basal glide between 50 and 
650K [13[ [I4j, and two alternative mechanisms with the change of temperature are 
proposed as the motion of a screw <a> dislocation on the prismatic plane through the 
cross slip process from the basal plane pL5j . According to this, in the low temperature 
region, the motion of screw dislocations is jerky and the rectilinear segment glides over a 
long distance before it is locked at the position of the basal plane. On the other hand, in 
the higher temperature region, the jog-pair between neighboring basal planes is formed 
and the glide is controlled by this mechanism. 

The core structures of screw and edge <a> dislocations and their gamma surface 
energies are studied by using the density functional method (DFT) and the EAM 
interatomic potentials to reveal properties of dissociation and the Peierls stresses at 
T = OK [16j. They found that the dissociated core structures predicted by the EAM 
potential m are in good agreement with the DFT results, although there are some 
discrepancies in the stacking fault energies(SFEs). Recently the core structure and the 
glide properties of the <a> dislocations have been studied in further detail by DFT [18]. 
and values of the Peierls stress for the basal and the prismatic glide have been predicted. 
The effects of the presence of solutes on the cross-slip softening of magnesium alloys are 
studied by modeling the double-kink nucleation (or the jog-pair) mechanism in the 
high temperature region m The DFT data on the interaction of prismatic and basal 
dislocation cores and the solutes are incorporated into the model and cross-slip stresses 
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are predicted as a function of solute concentration and temperature. 

In this paper, we track directly the cross-slip process of a screw <a> dislocation 
from the basal to the prismatic plane in pure magnesium using the DFT and the 
molecular dynamics calculations. 

2. Simulation methods 

In the DFT calculations, the electronic structure calculations and the structure 
relaxations by force minimizations are performed using the Vienna Ab-initio Simulation 
Package (VASP) [201 l2T1j with the projector augmented wave method and ultrasoft 
pseudopotentials. The exchange correlation energy is calculated by the generalized 
gradient approximation (GGA) with the Perdew-Burke-Ernzerhof function [22] . The 
Methfessel-Paxton smearing method with 0.2-eV width is used. The cutoff energy for the 
plane-wave basis set is 280 eV. Structural relaxation is terminated when the maximum 
force acting on the movable degrees of freedom becomes less than 10 rneV/A. 

Periodic quadrupolar array configuration of dislocation core shown in figured] is used 
to calculate core structures and their energies. The calculation cell is L[0110] x L[0002] x 
[2110] /3 with lattice constants ao = 3.1942 Aand Co = 5.1742 Awhich are derived from 
DFT calculation of a unit cell relaxation. These three orthogonal directions are referred 
to as x,y and z hereafter. Two system sizes of L = 8 and L = 12, consisting of 256 
and 576 atoms respectively, are used. In either case, the thickness of the cell is equal 
to the Burgers vector length b = | (a) |. A uniform shear strain is applied to the cell 
to cancel out the strain induced by the dipole moment of two cores [23]. For both x 
and y edges of the cell, b/2 is added to their z component. A 1 x 1 x 16 Monkhorst 
Pack k-point mesh is used for both sizes. Convergence of the energy difference between 
different core structures with respect to the increasing k-point mesh number and cutoff 
energy is confirmed. The numerical errors for the dislocation energy coming from the 
k-point mesh, the energy cutoff and the termination of the structural relaxation are 3, 
0.5 and 2 meV/6, respectively. The combined numerical error, excluding the finite size 
effect and the error coming from the GGA, is 4 rneV/fe. 

The core structure dissociated on the basal and prismatic plane is calculated by full 
relaxations. The initial configurations are obtained from the isotropic linear elasticity 
solutions for the displacement field of undissociated cores. Ideally one should use 
anisotropic linear elasticity solutions but the subtle difference between them does not 
affect the final relaxed structures, as long as periodic boundary conditions are used and 
no atoms are fixed to the positions assigned by the solution. To ensure that the core 
dissociates on the desired plane, we use linear combinations of two displacement fields, 
one is the original displacement field and the other is translated by the [0110] and [0002] 
direction for the basal and prismatic orientations, respectively. After the relaxation, 
the displacement field is calculated from the relaxed configuration and again the linear 
combination of the translated and untranslated displacement field is calculated and 
used as the initial configuration for the next DFT relaxation, to confirm that further 
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dissociation does not occur. In both basal and prismatic orientations, two core structures 
are obtained by initially placing the center of the core in the center or on the edge of a 
parallelogram, as shown in figure [2] Hereafter we refer to the four core structures of a 
screw dislocation shown in figure [2] as basal-0, basal-1, prism-0 and prism-1. 

The Peierls energies for the prismatic slip and cross slip are calculated using the 
drag method. For the prismatic slip, the initial and final state for the drag method 
is prism-0 and prism-1, respectively. For the cross-slip process, basal-0 and prism-1 
are used as initial and final configurations. In the calculations of the drag method, 
the intermediate configurations are firstly obtained by the linear combination of the 
initial and final configurations with a varying coefficient, and two Mg atoms which are 
closest to the center of the core in both initial and final configurations are fixed in the 
Burgers vector direction, then all the other degrees of freedom are relaxed using the 
DFT calculations. The choice of the two fixed atoms in each slip process is shown in 
figure |3j The drag method is easy to implement but in some cases fails to locate 
the saddle point m, and the failure is indicated by a discontinuous change of atom 
positions along the reaction coordinate. Therefore we have checked the atom positions 
and confirmed that they change continuously in all the cases investigated. 

Since the cross-slip process involves both the migration and the change of the plane 
of the dislocation dissociation, it is important to quantify both the positions and the 
structure of the dislocation core in the drag method calculations. Therefore, we invented 
a method described as follows. Let us denote the displacements in the Burgers vector 
direction of four atoms around the core normalized by the Burgers vector length by 
«i, u 2 , u 3 , n 4 as shown in figure [4j and U \ 2 — U\ — u 2 , and so on. There is one restriction 
among the four quantities, that is, U\ 2 + M 23 + W 34 + n 4 i = ±1. This leaves three degrees 
of freedom for the core structure. We define three quantities as follows: 

U x = m 34 - a 1 2 

Uy = U 23 — U41 

U s = W23 + u Al — m 34 W12 ( 1 ) 

These three quantities roughly correspond to the core position in the basal and prismatic 
direction, and the core structure, respectively. The positive and negative values of 
U s indicate dissociation on the prismatic and the basal planes, respectively. These 
quantities are used as reaction coordinates of the various slip processes. 

Molecular dynamics simulations are performed by the LAMMPS code [25] using 
the interatomic EAM potential mi- The time step is taken as lfs and the total number 
of time steps is at the most lxlO 6 . The simulation system is a Mg perfect crystal 
with the dimension of 130(x)xl22(r/)x637(,s)A( L x xL y xL z = 24x24x200). A screw <a> 
dislocation is introduced in the center of the system with the dislocation line in the 
z:[2110] direction. The AtomEye [26] is utilized for visualization of dislocations and 
defects in solids. The dissociated core structure of a dislocation in an hep crystal is 
well represented by the coordination number representation (figure [5]) in the same way 
as the case of fee crystals [2?]. When the temperature is above 150K, the visualization 
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becomes difficult due to fluctuations of atom positions. In the case of 300K and 400K, 
the ensemble average is taken from samples in the output interval of less than lps to 
visualize the dislocation structures.The periodic boundary condition is applied in the 
z direction and the solid side wall perpendicular to the x direction is displaced in the 
0 direction at the constant velocity of 0.01 — lOm/s to apply only the shear stress on 
the dislocation to promote a cross slip from the basal plane to the prismatic plane. In 
this system, a shear stress on the basal plane is inhibited. The spherical solid obstacle 
with the diameter of 6A is placed on the basal plane to study its effect on the cross¬ 
slip process. The obstacle is modeled as the aggregate of atoms on which no external 
forces are exerted. The dislocation that drifts on the basal plane due to the temperature 
cannot penetrate into the obstacle region. 

3. Results 

3.1. Results for the DFT calculations 

Figure |6] shows the core structures of a screw dislocation obtained by the DFT 
calculations of L = 12 case. The results, referred to as basal-0, basal-1, prism-0 and 
prism-1, are shown by the differentiated displacement map. Arrows are normalized in 
such a way that the differentiated displacement is 6/4 when the length is equal to that 
of the edge between atoms. The core structure of basal-0 agrees well with the previous 
DFT calculations [16]. The basal-1 is also stable, and the energy difference between the 
basal-0 and the basal-1 structure is less than 1 meV/b in both L = 8 and L = 12 cases. 
A precise estimation of the small barrier is out of the scope of this paper and it is suffice 
to confirm that the basal slip is extremely easy compared with other slips. 

As for the prismatic core structures, we found that the prism-0 and prism-1 
structure are both stable in the DFT relaxations. This stability is probably ensured 
by the exact mirror and/or inversion symmetry of the core which is strictly preserved 
in the DFT relaxations. The energy difference between the prisnr-0 and the prism-1 
structure is 0.5 rneV/ b and 1.5 meV/b for the L — 8 and L — 12 cases, respectively. 
The difference is smaller than the numerical accuracy and is expected to be sensitive to 
the system size and the boundary conditions. Thus again we only conclude that once 
the core structure changes from the basal to the prism-0 or prism-1, it easily glides on 
the prismatic plane. The Peierls stress for the basal and the prismatic glide have been 
calculated using DFT by Shin and Carter [18]. Their estimates are 0.6 and 35.4 MPa 
for the basal and the prismatic glide, respectively. 

Figure [7] shows the Peierls energy for the transition process from the basal-0 to 
the prism-0 core. Four atoms around the dislocation core are fixed in the Burgers 
vector direction during the intermediate states. Owing to the inversion symmetry, this 
is equivalent to fixing two atoms for the cases without the symmetry. For both L — 8 
and L — 12 cases, the stability of the prism-0 core is marginal. The maximum force 
acting on the atoms in the prism-0 core structure is less than 10 meV/A, and the 
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energy landscape is nearly flat around the prism-0 core structure. 

It is reported in Ref. [18] that the prism-0 core structure is actually unstable and 
transforms to the basal-0 structure. The stability of the prism core structures directly 
affects the temperature dependence of the cross-slip behavior; At the lower temperature 
region, experimental observations indicate that the dislocation cross-slips several atomic 
layers before reverting to the basal core structure. To explain this behavior, it has been 
assumed that the prism-to-basal core structure change requires some activation energy 
[15]. Since the results of the DFT calculations indicate that there is no activation energy 
for the prism-to-basal core structure change, we speculate that the low-temperature 
cross-slip behavior is a result of some inertia effect. 

The energy difference between the basal-0 and the prism-0 core structure is 47.4 
rneV /b and 45.4 rrieV /b for L = 8 and 12 cases, respectively. For comparison, figure [8] 
shows the corresponding Peierls energy calculated using the EAM potential for L — 12, 
24, 48 and 96. The energy of the prism-0 structure is used as the origin for clarity. One 
can see that the prism-0 core is stabilized in this case owing to the artificial minima in the 
stacking fault energy landscape in the prismatic plane, as reported in |28j. The Peierls 
barrier for the transition process from the basal-0 to the prism-0 is estimated as 16 
meV/5 from the L = 96 case, which is about one third of the DFT value. The finite size 
corrections in the EAM case is greater than in the DFT case, because the dissociation 
width of the basal core is much wider in the EAM case (about 12 A) compared with 
the DFT case (about 7.7 A), although the core shape parameter U s is almost the same. 
Figure [9] shows the Peierls energy for the cross slip calculated using DFT for the cases 
of L = 8 and 12. The initial and final configuration is basal-0 and prism-1 structure, 
respectively. The energy difference between the two structures is found to be 48.9 meV 
and 46.0 rneV/5 for the L = 8 and L = 12 cases, respectively. 

If a logarithmic interaction energy between dissociated partials and their mirror 
images is assumed, the finite size correction for the energy of dissociated core structure 
with the width d is proportional to d 2 / L 2 in the leading order. Its coefficient, including 
its sign, depends on the arrangement of periodic images of the core and the elastic 
anisotropy. Local strain induced by the I\ stacking fault of the width d also contributes 
to the finite size correction on the Peierls energy which is proportional to d 2 /L 2 . To 
estimate the Peierls barrier for the cross-slip more precisely, additional DFT calculations 
with higher precision and larger sizes are carried out. The energy difference between the 
basal-0 and the prism-1 core structure is calculated for L = 8, 12 and 16 cases, using 
the k-point mesh 2 x 2 x 32, 2 x 2 x 32 and 1 x 1 x 32, respectively. The structural 
relaxation is terminated when the maximum force is less than 2 meV/A. The overall 
numerical error is 2 meV/fe. Figure flOl shows the energy difference plotted against L~ 2 . 
The expected linear dependence is clear and we estimate the large L limit value as 
61.4 ± 2.0 meV/5. This is our final estimate for the Peierls barrier of the cross slip. 

In the cross-slip process, the core transforms from the basal to the prismatic 
structure and moves in the prismatic direction. It is an interesting question whether the 
processes of the transformation and the movement are simultaneous or sequential. If it 
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is sequential, the applied stress is of little help for the first transformation process. Thus, 
a high temperature or an extremely strong applied stress is required to initiate this first 
process. In order to see how the entire process proceeds, two dimensional Peierls energy 
landscape as a function of both U s and U y are estimated using the L — 12 case as shown 
in figure dU Beside the intermediate structures between the basal-0, the prism-0 and 
the prism-1 structures which are already presented, several additional core structures 
are investigated to cover the entire region. In these cases, the atomic configurations 
are prepared by the linear combinations of the basal-0, the prism-0 and the prism-1 
core structures, and the four atoms around the basal-0 core position are fixed in the z 
direction in the energy minimization. The gradient of the energy landscape is calculated 
from the forces acting on the fixed atoms and shown by arrows to help estimating the 
energy contour lines. The estimated minimum energy path from the basal-0 to the prism- 
1 structure, shown by the bold arrow, is close to a straight line, indicating that the core 
structure transformation and the movement of the core take place simultaneously and 
the process can be activated by a modest applied stress. The critical resolved shear 
stress r c can be roughly estimated by 

t c ~ AE/dbi 2 , (2) 

where A E is the Peierls barrier per b and d is the distance of the core positions between 
the initial and the final state. With A E = 61 meV/6, b = 3.19Aand d = 1.29A, the 
critical value is estimated to be about 740 MPa. 

For comparison, figure [12] shows the corresponding Peierls energy landscape 
calculated using the EAM potential for the L — 48 system containing a single 
dislocation. The boundary atoms are surrounded by vacuum region and are fixed in the 
x and y directions. The displacements in the z direction of these atoms are optimized for 
the prism-0 structure, and these displacements are used for all other cases. The barrier 
for the minimum energy path shown by the white arrow in £gure[[2]is 16 rneV/6. Figure 
[T3l shows the same energy landscape when a uniform shear strain of 0.01 is applied on 
the x plane in the z direction. The corresponding shear stress observed in each case 
is in a range between 207 MPa to 210 MPa. The minimum energy path moves toward 
the positive U y direction and the barrier reduces to 7 meV/6. ^Frorn these results, the 
critical resolved shear stress for the EAM potential is estimated to be about 400 MPa. 

Since the prismatic cross slip is a thermally activated process, it is crucial to 
estimate the activation enthalpy of the process in which a pair of jogs nucleate from the 
basal dislocation line towards the prismatic direction. The jog pair nucleation enthalpy 
can be estimated from a line tension model EMEU; 

Elt = J2 v ( 7 h) + \ K (jh - Vi- 1) 2 + b 2 hr]iT zx , (3) 

i 

where rji is a reaction coordinate of dislocation migration at the atomic layer i, V(rj) 
is the Peierls energy per b , K is a stiffness constant of the dislocation line, h is the 
dislocation migration distance and r zx is the shear stress. The three terms in the RHS 
of the equation [3] describes the energy contributions from the Peierls energy, the increase 
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of the length of the curved dislocation line and the work done by the movement of the 
dislocation, in their respective order. The harmonic string approximation described by 
the second term is valid unless some part of the dislocation line becomes edge-like, in 
which case the long-range interaction between dislocation segments becomes important. 
As we will see later, the dislocation line remains close to the screw direction even when 
a kink pair is generated and the harmonic string approximation is valid in the present 
case. 

We set p — 0 and p = 1/2 for the basal-0 and the prism-1 core structures, 
respectively. We consider the jog-pair nucleation in the high temperature region, 
in which the prismatic dislocation core reverts to the basal core at the next basal 
plane. Thus p — 1 corresponds to the basal-0 core on the adjacent plane, and 
h = Cq/2 = 2.587A. The coefficient K is estimated by the DFT calculation of two-layer 
system which is composed of two L — 8 configurations with p\ =8/16 and p 2 = 7/16 
|29j . Our estimate is K = 8.15 ± 1.03eV, or in terms of the line tension, 5.5 to 7.0 nN. 
The DFT data for V(p) is fitted by a function Xp=i C/(l — cos(2jnp))/2 with C\ = 55.95 
meV, C 2 = 19.89 meV, C 3 = 5.20meV and C 4 = 1.28 meV, using the values in figure [9] 
scaled by a factor to reproduce the final estimate of the barrier 61.4 meV/5. 

The activation enthalpy is estimated by the minimization of the equation [3] to be 
1.4 to 1.7 eV for the r zx = 0 MPa case and 1.1 to 1.3 eV for t zx = 100 MPa case, 
considering the numerical error for the Peierls potential and K. Figure fl4l shows the 
saddle point configuration of the jog-pair nucleation for the cases of r zx = 0 MPa and 
100 MPa. The swept area, which is equal to the derivative of the activation enthalpy 
with respect to the applied stress divided by b, is about 12 b 2 and 8 b 2 for the 0 MPa and 
100 MPa cases, respectively. The experimentally observed activation enthalpy for the 
prismatic slip is reported to be 1.0 to 1.4 eV between 315 and 664K, and 1.63eV at 590K 
|13j . Thus our result is in reasonable agreement with these experimental observations. 
The activation area observed in experiments is 6 to 12 b 2 at 300K [15], which is also in 
agreement of our results. 

3.2. Results for the molecular dynamics method 

The core of a screw <a> dislocation on the basal plane is dissociated into two partial 
dislocations as shown in figure [5] With the present system size, the same structure 
is obtained by equilibrating with the system temperature of T = 0.01K regardless of 
the boundary conditions. The widths of the stacking fault between these two partials 
are 12.38 and 12.4lA at T = 0.01 and 10K, respectively. Since an <a> dislocation 
on the basal plane is easily moved by the disturbance of temperature, the dissociated 
core structure cannot be determined precisely above 50K due to the meanders of the 
dislocation line. 

Figure [15] shows the temperature and strain rate dependence of the shear stress 
vs strain curves for the cross slip of a screw <a> dislocation from the basal to the 
prismatic plane. The shear zx-component of stress is the averaged value over the entire 
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simulation region. Firstly, the effect of strain rate is studied on the critical shear stress 
for the constriction and the onset of a cross slip. The strain rate is varied from 2.5xl0 6 
to 2.5x10 9 [1/s] at T = 0.01K. It is seen from figure H5l that the critical shear stress 
is largely affected by the strain rate and is converged at 2.5x10'[1/s]. There is still 
a large difference from the experimental value of less than nearly 1 [1/s], but this 
value of strain rate is in the relatively lower region with the present-day molecular 
dynamics simulation. Therefore, the converged strain rate of 2.5x10'[1/s] is used to 
study the temperature dependence. From figure [HI the critical shear stress decreases 
with increasing temperature to 50K, where the critical strain also decreases. Above 50K 
to 400K, the critical shear stress continues to decrease with increasing critical strain. 
This indicates that the critical shear stress for the cross slip decreases with increasing 
temperature, while the time to the constriction increases. In other words, the material 
softens with increasing temperature above 50K. In order to exclude this effect, the 
critical shear stress, which is scaled by the slope of the stress vs strain curve, that is the 
shear stress fi, corresponding to each temperature, is shown in figure flfil as a function of 
temperature. It is seen that the scaled critical shear stress reduces down to 50K due to 
the thermal activation effect. The scaled critical shear stress, or the critical shear strain, 
increases above 100K, leading to the delay of the cross slip. The reason for the delay of 
the cross slip or the constriction is that the fluctuation motions of the dislocation line on 
the basal plane become large as the temperature increases, and this motion leads to the 
delay of the stress build-up on the dislocation line. The critical shear stress is expected 
to behave r c = r*(T) — AT , where r*(T) denotes the critical shear stress at which the 
activation free energy for the cross slip becomes zero, and A is some positive constant. 
Above 50K, t*(T) is expected to increase rapidly, leading to the overall increase in r c . 

Figure [171 shows the cross-slip sequence of a screw dislocation after the constriction 
at T = 0.01K. It is seen from the figure fTTlfa) that the constriction occurs at two 
places since the dislocation is almost motionless on the basal plane at this temperature. 
Although the critical shear stress for constriction is high, the dislocation moves easily 
on the prismatic plane once it is constricted. In figure [IW b). the cross-slipped part of a 
dislocation consists of rectilinear segments with a rather long linear segment in the front 
center gliding over a long distance. This corresponds to the locking-unlocking process 
as proposed in [15] . A close look at the bowed-out segment of a dislocation line in figure 
[171( b) indicates that the curved part consists of slightly dissociated linear dislocation 
lines with connecting jogs between neighboring basal planes. This suggests that the 
jog-pair mechanism |15L [TO] , which applies to the case in the high temperature region, 
works locally even at this low temperature along with the jerky motion. 

Figure [18] shows the glide of the cross-slipped dislocation on the prismatic plane 
at different temperatures. All the figures are the snapshots taken at 5ps after the 
constriction of the dissociated screw dislocation on the basal plane. It is seen that 
the glide motion of the dislocation is apparently retarded as the temperature rises 
above T = 50K, where the curved bowed-out dislocation line becomes smooth. This 
means that the length of the linear segments consisting the curved line becomes smaller. 
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Also, at finite temperatures, plateau steps appear on the curved part of the bowed-out 
dislocation line. These plateau steps appear even at T = 10K where the center of 
the bowed-out part travels a long distance. These plateau segments are dissociated on 
the basal plane and the motion of a dislocation on the prismatic plane is retarded by 
these segments. In particular, the development of this plateau segment is noticeable at 
T = 400K as seen in figure fTSf f h in which the linear part is fully dissociated on the basal 
plane. Because of the high strain rate of the molecular dynamics method, the formation 
of this dissociated linear plateau segment observed at T = 400K may actually be shifted 
to the lower temperature region and corresponds to the transition of the motion of a 
screw dislocation around T = 250K observed in the experiment |15j. Figure fl9l shows 
the temperature dependence of the area swept by the bowed-out dislocation on the 
prismatic plane after 5ps from the constriction. The reduction of the swept area is 
noticeable above T = 50K by the formation of the plateau steps in the bowed-out part 
of the dislocation line. The constricted part of the dislocation enlarges due to the 
increase of temperature above T = 50/1, in particular, around T = 300/1. The dip 
in the swept area at T = 400/1 is caused by the large dissociation of the dislocation 
at the plateau step as observed in figure U8lf). It should be emphasized that this 
result indicates only the early stage of the cross-slip process and does not mean the 
glide property of an <a> dislocation on the prismatic plane. 

Figure [20] shows the shear stress vs strain curve when the hard sphere obstacle with 
the diameter of 6A is placed on the same basal plane in which the screw dislocation is 
set. When the dislocation is contacted with the obstacle, the local stress is concentrated 
at the contacting part and a large reduction of the critical shear stress for the cross slip 
is observed at T = 50 and 150K. The effect of the strain rate becomes large compared 
with the no obstacle case where a large reduction of the critical shear stress is obtained 
as the strain rate decreases. The critical shear stress converges at the present shear stress 
of 2.5x10'[1/s] as in the case of no obstacle. At higher temperatures, near the room 
temperature, the reduction of the critical shear stress is erratic because the dislocation 
line meanders on the basal plane. A large reduction of stress is not always obtained as 
in the case of low temperatures. In figure l2TTfa) at 50K, it is seen that the cross-slipped 
dislocation is pinned by the obstacle and the motion on the prismatic plane is retarded. 
At higher temperatures, the dislocation is found to be easily unpinned by the obstacle. 

4. Discussion and conclusions 

The Peierls energy and the activation enthalpy for the prismatic cross slip of a screw 
<a> dislocation in Mg are predicted using the DFT for the first time. The energy 
difference between the basal and the prismatic core, which is the Peierls barrier for the 
prismatic cross slip, is 61.4 ± 2 meV per Burgers vector length. On the other hand, 
in the widely used EAM potential, the Peierls potential has a peak at the intermediate 
structure between the basal and the prismatic core, and the barrier height is 16 meV 
per Burgers vector length. 
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The activation enthalpy and the activation area of the jog-pair mechanism, which 
operate in the high temperature range, are estimated using the line tension model based 
on the DFT results. The activation enthalpy and the activation area are estimated as 
1.4 to 1.7 eV and 8 h 2 to 12 h 2 , respectively. They are in reasonable agreement with 
experimentally observed values of 1.0 to 1.6 eV and 6 b 2 to 126 2 , respectively [13] H5j . 

The Peierls stress estimated from the Peierls barrier is about 740 MPa and 460 
MPa for the DFT and the EAM cases. However, it is known for the case of body centered 
cubic metals that the atomistic prediction of the Peierls stress becomes far greater than 
the experimentally observed yield stress in the low temperature region, and the cause 
of the discrepancy has not yet been fully identified [32] . In pure Mg, the reported yield 
stress for the prismatic slip in the low temperature region is on the order of lOOMPa. 
Thus the discrepancy is also present in the Mg case. 

The estimate of the activation enthalpy and the relevant core structures are crucial 
to predict the alloying elements which catalyze cross slips. By comparing the solution 
energies of the solute atom near the basal and the prismatic core, one can calculate 
the reduction of the activation enthalpy by the solute atom. If it is comparable to the 
original activation enthalpy of ~ 1.5 eV, the prismatic cross slip is expected to occur 
at the lower temperature. The results of the MD simulations indicate that the thermal 
fluctuations of the basal core hinder the activation of the prismatic cross slip, and an 
obstacle for the basal slip suppresses the fluctuation and reduces the CRSS for the 
cross slip. This result indicates that a solute atom which has either a repulsive or an 
attractive interaction with the basal core also reduces the activation enthalpy of the cross 
slip. Further DFT studies on the solute-dislocation interactions and the development of 
interatomic potentials for Mg accounting for the proper dislocation properties such as 
the Peierls barrier are expected. 
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[2110]/3 

Figure 1. Rectangular cell used in the DFT calculation, seen from the Burgers 
vector direction. The white circles are Mg atoms and the two signs ”+” and mark 
the positions of two screw dislocations with opposite helicities, which form a periodic 
quadrupolar array. 



Figure 2. Four core structures of a screw dislocation investigated using DFT. 
The ellipses show the position of the core, either on the edge or at the center of the 
parallelograms, and the planes of dissociation. 
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Figure 3. 

black disks, 
slip process. 


Two Mg atoms fixed in the calculations of the drag method are shown by 
FIX Red and blue ellipses are the initial and final configuration of each 



Figure 4. The four differentiated displacement Uij around a dislocation core used to 
quantify the core position and core shape. See the main text for more details. 
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Figure 5. A dissociated screw <a> dislocation on the basal plane in the coordination 
number representation. 
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Figure 6. The four core structures of a screw dislocation on the basal and prismatic 
planes obtained by the DFT calculations. White circles show the in-plane positions of 
each atom, and arrows show differentiated displacement map. The black squares mark 
the center position of the core, and the dotted lines are the stacking faults where the 
differentiated displacement is greater than 6/2. 
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Figure 7. Peierls energy of the transition process from the basal-0 to the prism-0 
core structure calculated by the DFT. 



Figure 8. Peierls energy of the transition process from the basal-0 to the prism-0 
core structure calculated using the EAM potential. 



Figure 9. Peierls energy of the cross-slip process from the basal-0 to the prism- 1 
core structure, calculated by the DFT. 
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Figure 10. Energy difference between the basal-0 and the prism- 1 core structure 
plotted against L~ 2 . The arrow indicates the large L limit estimated from a linear fit. 



Figure 11. Two dimensional Peierls energy map as a function of U s and U y , estimated 
from the DFT results. The number above each square shows the Peierls energy (meV /b) 
calculated by the DFT. Gradient of the energy landscape is calculated from the forces 
acting on the fixed atoms and shown by the arrows to help estimating the contour 
lines. The bold arrow indicates the minimum energy path from the basal-0 core to the 
prism-1 core structure. 
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Figure 12. Two dimensional Peierls energy landscape as a function of U s and U y , 
calculated using the EAM potential for the L = 48. The bold arrow indicates the 
minimum energy path from the basal-0 core to the prism-1 core structure. 
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Figure 13. Two dimensional Peierls energy landscape when 1% shear strain is 
applied on the x plane in the z direction, calculated using the EAM potential for the 
L = 48 case. The bold arrow indicates the minimum energy path from the basal-0 
core to the prism-1 core structure. 
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Figure 14. The saddle point configurations for the jog-pair nucleation, derived from 
the line tension model for the case of applied shear stress OMPa and lOOMPa. 
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Figure 15. Shear stress vs strain curves for the cross-slip process at different strain 
rates and temperatures. 
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Figure 16. Temperature dependence of the scaled critical shear stress for the 
cross-slip process. The critical shear stress is scaled by the shear modulus /i for each 
temperature. 



(a) lps after the constriction. 



(b) 3ps after the constriction. 



(c) 5 ps after the constriction. 


Figure 17. Glide of a screw <a> dislocation on the prismatic plane at T = 0.01K, 
where the side view (left) and the plan view (right) are shown. 
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(c) T = 100K 


(d) T = 150K 



(e) T = 300K 


(f) T = 400K 


Figure 18. Glide of a screw <a> dislocation on the prismatic plane at different 
temperatures. The side views are taken at 5ps after the constriction. 
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Figure 19. Temperature dependence of the area swept by the cross-slipped 
dislocation on the prismatic plane after 5 ps from the constriction. 
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Figure 20. Reduction of the critical shear stress by the solid sphere obstacle on the 
basal plane. The diameter of the obstacle is 6A. 



(a) T = 50K 



(b) T = 150K 


Figure 21. Glide of a cross-slipped screw <a> dislocation on the prismatic plane 
when the obstacle is placed on the basal plane. 























